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For continuous weak measurement of qubits, we obtain exact expressions for weak values (WVs) 
from the post-selection restricted average of measurement outputs, by using both the quantum- 
trajectory-equation (QTE) and quantum Bayesian approach. The former is applicable to short-time 
weak measurement, while the latter can relax the measurement strength to finite. We find that 
even in the “very” weak limit the result can be essentially different from the one originally proposed 
by Aharonov, Albert and Vaidman (AAV), in a sense that our result incorporates non-perturbative 
correction which could be important when the AAV’s WV is large. Within the Bayesian framework, 
we obtain also elegant expressions for finite measurement strength and find that the amplifier’s noise 
in quantum measurement has no effect on the WVs. In particular, we obtain very useful result for 
homodyne measurement in circuit-QED system, which allows for measuring the real and imaginary 
parts of the AAV’s WV by simply tuning the phase of the local oscillator. This advantage can be 
exploited as efficient state-tomography technique. 

PACS numbers: 03.65.Ta,03.65.Yz, 42.50.Lc,42.50.Pq 


I. INTRODUCTION 


The concept of weak value (WV) was introduced by 
Aharonov, Albert and Vaidman (AAV) in the seminal 
paper jlj and was elaborated with extended discussions 
H D. The most interesting and surprising feature of the 
WV is that it can exceed the range of eigenvalues of the 
observable. Obviously, the appearance of this nonclassi- 
cal or strange WV violates our common knowledge that 
the quantum average must be bounded by the extremum 
of eigenvalues. Therefore its interpretation has been a 
subject of confusion §!• Despite of this, the strange 
WV has been experimentally confirmed [|] |). 

In addition to theoretical curiosity, the concept of WV 
is very useful to facilitate to resolve and generate quan¬ 
tum paradoxes, such as the Hardy’s paradox §@,the 
superluminal transport in optical fiber 0 or in vacuum 
|Q , the momentum-disturbance relationship [jld) and lo¬ 
cally averaged momentum streamlines 0® in two-slit 
interferometer, and the more recently proposed quantum 
Cheshire Cat ||l6| |l8|. Even more, the WV has novel ap¬ 
plications in techno ogical aspects, e.g., for weak signal 
amplification and sensitive estimation of unknown small 
parameters |[l9|, fiol , and for efficient quantum state to¬ 
mography |21] |2.f . 

In AAV’s original work, the WV of spin-i particle was 
proposed as (a z ) w = Re Jd* = Re ( cr ^)’ where 1 i>i) 
and \ipf) are the pre- and post-selected states (PPS). 
Note that the key quantity (we term it in this work as 
AAV’s WV), <7* = Jq i is a complex object. Its 
interpretation is thus subtle and has caused debates ]2] 
|f|. Simpler understanding for WV is from the practical 
measurement point of view j?|, || , which is to be termed 
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as measured WV or simply as WV in this work. The mea¬ 
sured WV corresponds to a PPS restricted average, say, 
a sub-ensemble-average of the weak measurement results, 
with partial data survived from the post-selection. 

In most of the existing WV studies, the weak measure¬ 
ment is implemented via a unitary interaction model by 
two coupled subsystems, i.e., the measured system and 
the measuring device, while the device is described by a 
quantum Gaussian probe wavefunction with continuous 
variables which is subject to extra projective measure¬ 
ment for the “pointer” coordinate or momentum. For 
this type of model and using a conditional-probability 
formulation the WV can be re-derived as Q : 

(a z ) w = = J ^-P x (f ; i)dx 

= Re(cr^) + 0(e 2 ), (1) 

where e is the unitary interaction strength which deter¬ 
mines also the average shift of the pointer caused by the 
system eigenstate. P x (f;i) is the success probability of 
post-selection with |"0/}, following the x-conditioned sys¬ 
tem state after the weak measurement. Note that this 
post-selection conditioned average is also the method of 
data analysis in experiments Jl], |^| . 

One may note that, differing from the unitary inter¬ 
action model mentioned above for weak measurements, 
there exists alternative class of models and real systems, 
which has been extensively studied for qubit measure¬ 
ment and control, especially in solid-state setup using 
such as the quantum point contact (QPC) (25| |27|, the 
single electron transistor (SET) |^8[ ^9|, or the homodyne 
detection in circuit quantum-electrodynamics (cQED) 
f{()| fi~if. However, despite particular importance, this 
type of weak measurement has not yet received much 
attention in the context of WV studies. Its main dif¬ 
ference from the unitary interaction model is that the 
“device” in this type of measurement scheme is described 
“classically”, without involving the usual quantum wave- 
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function description. Instead, this type of weak mea¬ 
surement is characterized by quantum trajectories — the 
stochastic evolutions of the system state conditioned on 
the (continuous) measurement outputs, which can be de¬ 
scribed by quantum-trajectory-equation (QTE) or quan¬ 
tum Bayesian rule. In Ref. |35|, applying QTE and using 
the conditional-probability method, the canonical AAV 
result of Eq. ([I]) was recovered for continuous weak mea¬ 
surement of qubit. 

In present work, along the line of Ref. [^5|, we recalcu¬ 
late the WVs for this type of weak measurement. Beyond 
Eq. (@), we obtain the following better result: 

f(x)i Re(gg,) 

e l + 5(kS,| 2 -l)- U 

Similar as in Eq. ( 0 ), the scaling parameter e in this ex¬ 
pression is the shift of the average output current caused 
by the basis state of the qubit. Notably, in the denomina¬ 
tor, g characterizes the measurement strength and takes 
a simple form g = 7 dt with 7 the measurement rate and 
dt the measurement time. 

Remarkable modification of Eq. ( 0 ) over Eq. ( 0 ) is the 
“partial”-summation-type term in the denominator. Ac¬ 
tually, as long as the expansion e~ 2g ~ 1 — 2 g holds to be 
valid, Eq. ( 0 ) is the exact result of the WV. This conclu¬ 
sion can be proved by calculating f(x)i via the quantum 
Bayesian rule, which does not restrict the measurement 
strength to weak limit. We will see that the exact WV 
result for arbitrary measurement strength corresponds to 
a simple replacement g —> Q = (1 — e~ 2g )/2 in Eq. ( 0 ), 
while the measurement strength reads g = 7 t m , with t m 
now a finite time of measurement. 

The non-perturbative correction in Eq. cannot be 
neglected when < 7 ^ is large and should be taken into ac¬ 
count in the various applications of WVs. For instance, 
in the context of WV amplifications, it has been pointed 
out that some “nonlinear” correction to the linear result 
of AAV WV would set an upper bound for the amplifica¬ 
tion coefficient, based on similar results obtained within 
the framework of the unitary interaction model of weak 
measurements j|7] [ 50 ] . As a first and simple remark, our 
present work differs from those in that we consider dif¬ 
ferent type of weak measurement, which allows for em¬ 
ploying different formulations to obtain non-perturbative 
and exact WV expressions for the measurement of qubits 
under consideration, in particular obtaining useful result 
for the solid-state superconducting circuit-QED setup. 

The remaining part of this work is organized as follows. 
In Sec. II, we complete all the formal aspects of calcu¬ 
lating the WVs for the continuous weak measurement of 
solid-state qubits. In Sec. Ill we specify the study to 
solid-state circuit-QED setup to carry out useful result 
which allows for direct measurement of the real and imag¬ 
inary parts of the AAV WV, by simply tuning the local 
oscillator’s phase in the homodyne measurement. In Sec. 
IV we summarize the work with brief remarks. 


II. WEAK VALUES OF QUBIT 
MEASUREMENTS 


A. Continuous Measurement of Qubit 


There are several experimentally accessible systems 
for continuous weak measurement of qubits, includ¬ 
ing for instance the solid-state charge qubit measured 
by quantum-point-contact or single-electron-transistor 

g P9| , and the superconducting qubit in circuit-QED 
, [3 If measured via the cavity-field homodyne detec¬ 
tion [B0 B4[. Under reasonable simplifications, all these 
setups can be commonly described as follows. First, 
the qubit Hamiltonian is reduced to a two-state model 
H q = y a z , where a z = 11)(11 — |2)(2| and A = Ei — E 2 . 
Second, the measurement principle is based on the dis¬ 
tinct stationary current in the detector, Tip, which is 
associated with the qubit state |1) or |2). Moreover, 
of particular interest is continuous weak measurement 
for a superposition state \ipi) = a|l) + /3|2). In this 
case, the output current can be expressed in general as 
/(f) = I 0 + %^(cr 2 ) + £/(f), where I 0 = (h + h)/2 and 
A/ = |/i — I 2 \. The last term obeys the statistics of white 
noise, with ensemble average E[£/(t)£j(t')] = ^-5(t — t'), 
where Sq is the spectral density of the output current. In 
this work, we would like to adopt a reduced description, 
by denoting the measurement output as J(t) = -^ 7 ==. 

Then we have J(t) = 2^y(a z ) + £(t), where 7 = 
and E[£(t)£(t')] = 6(t — t'). 


1. QTE Approach 

Conditioned on the continuous outcomes of measure¬ 
ment given above, the qubit state evolution (in the rotat¬ 
ing frame with respect to H q ) is governed by the following 
Ito-type quantum trajectory equation (QTE) |uj 

dp = ^V[(j z ]pdt + y/~f'H[a z }pdW(t ). (3) 

The super-operators are defined as T>[A]p = ApA' — 
T{A^ A,p}, and TL[A\p = Ap+pA^ — TA[(A-\-A^)p]p. The 
Wiener increment is formally defined by dW(t) = £(f) dt, 
in experiment which should be extracted from the output 
current. Then, after a short-time (dt) measurement, the 
state is updated as p(t + dt) = p(t) + dp(t), in the qubit 
basis which reads 


P11 = P11 + 2^7(1 - W z ))Pn dW , (4a) 

P22 = P22 - 2 ^ 7(1 + (& Z ))P22 dW , (4b) 

P12 = P12 - 2 / 512(7 dt + ^/l((T z ) dW) , (4c) 

where ( a z ) = Tr [a z p(t)\ = pn — p 22 - Below we com¬ 
pare this result with the one from the quantum Bayesian 
approach. 
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2. Quantum-Bayesian Approach 


The quantum Bayesian approach, originally proposed 
by Korotkov |} 6 |, is based on the well-known Bayes for¬ 
mula in Probability Theory together with a quantum pu¬ 
rity consideration. The former is utilized to determine 
the diagonal elements while the latter is for determina¬ 
tion of the off-diagonal ones. In quite compact form, the 
quantum Bayesian approach updates the qubit state from 
p to p according to the following rule: 


Pu = pnPi(x)/AT{x), p 2 2 = P 22 P 2 {x)/N{x), 

P 12 = P\ 2 \JPl(x)P 2 {x)/M(x). (5) 

In these formulas, associated with the qubit state | j) 
(j = 1,2), Pj(x ) is the priori probability distribution 
“knowledge” for the stochastic outputs “x”. J\f(x) = 
pnPi(x) + p 22 P 2 {x) is a normalization factor. Note that 
the last rule (for off-diagonal element) is from a purity 
consideration. About the distribution “knowledge”, in 
most cases, the integrated output (denoted here by the 
stochastic variable “x”) satisfies a Gaussian statistics: 


Pj(x) 


, - exp 

v / 2 VD 



( 6 ) 


where Xi( 2 ) = ±e denotes the average output correspond¬ 
ing to qubit state 11 (2)), and D is the distribution vari¬ 
ance. 

As a more special case, let us consider a short-time 
measurement over dt. The measurement result (stochas¬ 
tic “x”) simply reads x = J(t)dt = 2^/ 7 y(a z ) dt + dW(t). 
From this, we can identify e = 2 ^/ 7 dt and D = dt. View¬ 
ing the small parameters e, dt and x, we first expand the 
exponential functions in the Bayesian formulas to obtain 


We may find that these two equations differ from the 
QTE (||) or (Q). Indeed, using the Bayesian rule in this 
manner, we obtain the Stratonovich-type QTE, but not 
the ltd- type To obtain the QTE of the Ito-type, one 
needs to add a term (J 7 / 2) (dP/dpij ) to each pa, where J 7 
is the factor before £(t) in each equation [g7|. Applying 
this conversion rule, one obtains 

pu = -P 22 = V 7 PiiP 22 .£(t), ( 10 a) 

P 12 = - 27 P 12 - 2 y / 7 y(a z )p 12 ^{t). ( 10 b) 

This result is the same of Eq. !) or (I- 

Below we make more comparisons. Inserting the above 
Bayesian expansion into the Ito calculus pij (t) = [pij (t + 
dt) — Pij (t)\/dt. one can obtain the Ito-type QTE under 
the Milstein algorithm f43| , for which we exemplify 

dpu = 4 V 7 PnP 22 dW - 87 / 211 ^ 22(^11 - P22) 

x {(dW) 2 - dt] . (11) 

Note that this result differs from dpu = 4 Wypnp 22 dIK, 
which corresponds to the so-called Euler algorithm |43|| . 

Alternatively, if we only keep dt and dW as small pa¬ 
rameters in the Bayesian expansions and use the calcu¬ 
lus Pij(t) = [pij(t + dt) — pij(t)]/dt, we find that expan¬ 
sion up to the first order of dt and dW corresponds to 
the Stratonovich-type QTE. However, if keeping terms in 
the expansion up to (dW) 2 , one arrives to the Ito-type 
QTE under the Milstein algorithm. This observation tells 
us that, in order to make the result from the quantum 
Bayesian rule be fully identical to that given by the Ito- 
type QTE, one should first make the expansion of the 
Baysian formula up to (dW) 2 , then set (dW) 2 = dt. 


P 11 = Pu 


P 22 — P 22 


0 

(> 




ex 


1 + D {r 


ex 


u-ir 


(7a) 


B. Weak Values 

1. QTE Approach 


(7b) The post-selection restricted average can be calculated 

using the following conditional-probability method: 


p 12 = P 12 


1 + !(—)* 

4 V 2 ZV 




-1 


(7c) 


/ dxxPi(x)P x (f) 
f{X)i ~ JdxPi(x)P x (f) 


( 12 ) 


Performing further expansion up to the second-order 
of e and x, one can prove their equivalence to Eq. (||) 
or ([|) as follows. Actually, using the Bayesian expan¬ 
sion and applying the calculus pij(t) = [pij(t + 4f) ~ 
Pij(t— ^)]/dt, we obtain for instance pu(t) = —f> 22 {t) = 
(2ex/D)pii(t)p 2 2(t), which gives 

PU = P22 = 87(<T~)pii/222 + 4 - v /7/9iiP22£(£) ■ (8 ) 

Similarly, for the off-diagonal element, we have 

P 12 = -47(cr 2 ) 2 pi2 - 2^/7(<7 2 )pi 2 £(£). (9) 


where x stands for the stochastic outcome of the weak 
measurement, with probability Pi(x) determined by the 
qubit state |V’i). Assuming | tpi) = a\l) + /3|2), we have 
Pi(x) = |a| 2 Pi(x) + |/ 9 | 2 P 2 (a)), where P\( 2 ){x) is the priori 
probability corresponding qubit state 1 1 ( 2 )). P x (f) is the 
success probability of the post-selection by state \4’f), fol¬ 
lowing the measurement outcome x. As explained above 
for the continuous measurement of qubit, the output cur¬ 
rent can be expressed as J(t) = 2^/7 (a z ) + £(t). For 
time interval dt , the output charge (integrated current) 
is x = J(t)dt , and the qubit state is updated via Eq. (||) or 
([Ij) by inserting dW = x — 2^y(a z )dt. Then, the success 
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probability of post-selection is P x {f) = {^f\p(t + dt)\ipf). 
More explicitly, we have 

P x (/) = P{1 + 7^(K| 2 -1) 

+ 2V7[Re«) - ( a z )]dW} , (13) 

where P = \{^f\^i)\ 2 and cr 2 = (V’/k*|V’i)/(V’y#i}, as 
defined in Sec. I (Introduction). Straightforwardly, we 
evaluate the numerator of Eq. (|l^) as 

Mi = P{[1 + ^dt{a z ){Re(cr z w ) - (<j 2 ))](|a| 2 xi + |/3| 2 x 2 ) 
+2^7(Re«) - {cr z ))[\a\ 2 (D + x 2 ) + \/3\ 2 (D + x 2 )]} , 

(14) 

and the denominator as 

M 2 = P[1 + ”tdt{\(T z w \ 2 - 1) - 47(cr 2 )df(Re(cr 2 ) - (cr 2 )) 
+2^7(Re(a 2 ) - {a z )){\a\ 2 Xi + |/3| 2 x 2 )]. (15) 

Accounting for |a| 2 + |/3| 2 = 1, |a| 2 — |/3| 2 = (cr 2 ) and 
xp 2 = ±2^/7 .dt = ±e, the WV takes then a very compact 
form as 

f {x)j _ Re(g~) 
e 1 + P(k 2 | 2 -1) ' 1 j 

Here, we have denoted g = 7 dt , which characterizes well 
the measurement strength. This result differs from the 
canonical AAV WV, f(x)i/e = Re(cr 2 ), by a “partial” - 
summation-type correction in the denominator. How¬ 
ever, as we will see in the follow ing , as long as the ex¬ 
pansion e~ 2g ~ 1 — 2 g holds, Eq. (|16|) is non-perturbative 
and exact. As shown in Fig. 1, even in the weak mea¬ 
surement limit (short time dt or small p), the role of this 
correction could be essential when a^ is large. This im¬ 
plies that in the various WV applications such as weak 
signal amplification, the non-perturbative correction in 
the denominator of Eq. © should be taken into account 

mm- 


2. Bayesian Approach 

Applying the quantum Bayesian rule we first calculate 
the post-selection probability, P x (f) = Tr [pfp(x)\, where 
Pf = is the post-selection state and p(x) the 

updated state from Eq. (||). The result simply reads 

Px{f) = [PfllPllPl{x) + P/ 22 /0 22 .P 2 (x)] 

+2Re(p} 12 p 12 ) V / Pi(z)P 2 0r). (17) 

Pij are the elements of the density matrix of the initial 
state p = which gives also the distribution func¬ 

tion Pi(x) in Eq. (|l2|) as Pi{x) = pi\P\{x) + p 22 P 2 (x). 
Straightforwardly, after completing a number of Gaus¬ 
sian integrals, the WV defined by Eq. (|l2]) is obtained as 
f(x)i = Mi/M 2 , with Mi and M 2 given by 

Mi = {pfllPuXi + Pf 22 p 22 X 2 ) + (xi + x 2 ) 
x Re(pj 12 p 12 ) exp [-(Si - x 2 ) 2 /( 8 D )] , (18a) 



Qln 

FIG. 1: (color online) Non-perturbative correction of Eq. ( 0 ) 
to the canonical AAV WV (blue line). In this plot, we specify 
the initial state as \ipi) = -^j(|l} + |2)) and the post-selection 

state as 11/7) = [(cos § + sin §)|1) + (cos § — sin |)|2)]/\/2, 
which result in the AAV WV crj = tan(0/2). We wee that, 
in contrast to the AAV WV, a different trend appears in the 
large crj regime. 


M 2 — ( PfiiPn + p/ 22 p 22 ) + 2 Re(p)- 12 pi 2 ) 

x exp [-(5i - x 2 ) 2 /(8 D)] . (18b) 


Since the Bayesian rule works for finite strength measure¬ 
ment, this result is then valid in general. 

To be more specific, let us consider an initial state 
| tjji) = ^7f(|l) + |2)) and parameterize the post-selection 
state | ipf) = di|l) + d 2 |2) with polar angle 9 as d\ = 
(cos | + sin |)/\/2 and d 2 = (cos | — sin |)/-\/2. Substi¬ 
tuting the two states into Eq. (|l8|), we obtain 


/(*)» 


Xi + x 2 
2 


sin 9 / x\ — x 2 \ 

1 + G cos 9 \ 2 J 


(19) 


where G = exp[— {x\ — x 2 ) 2 /(8 D)\. Similar as Eq. (]T|), 
this result is valid for arbitrary strength of measurement. 

Now consider the limit of weak measurement based on 
this result. If we take first the limit G —> 1, we find 
that the WV can diverge when 9 —> n, like the AAV 
WV, in present case which reads <r 2 = tan(0/2). This 
limiting order corresponds also to taking first the limit 
g = 7 dt — > 0 in the denominator of Eq. (pH) , then altering 
the post-selection to get divergent WV. However, for any 
real measurement, the measurement strength must be 
nonzero, which results in f(x)i/e behaving as shown in 
Fig. 1. The turnover behavior implies an existence of 
optimal post-selection state which maximizes the WV, 
and vanishing WV when the AAV WV diverges. The 
latter corresponds to G ^ 1 and 6 —> n in Eq. (p9|). 

Moreover, it is desirable to recast the general WV re¬ 
sult Eq. (jl^) into more compact form, like Eq. ( pb| ) in the 
weak measurement limit. In the Bayesian result, let us 
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identify 27,2 = ±2^7 t m = ±e and D = t m , where t m 
is the finite measurement time. Substituting these into 
Eq. (|l^) , after some algebra we get 

f(x)j _ Re(Q 

e 1+6(KP-1)' 1 j 

Corresponding to g = 7 dt in Eq. d!): the counterpart Q 
in this result has been generalized to Q = (1 — e~ 2g )/2 
with g = "ft m . Obviously, both results are identical for 
short time limit. This proves also that the QTE result 
of WV is exact in the weak measurement limit, having 
not lost anything in the “partial”-summation-type form 
of Eq. (@. 


3. Noisy Amplification 


For many realistic measurements such as in circuit- 
QED (to be addressed in Sec. Ill), the weak output sig¬ 
nal x must be amplified properly. As a simplified de¬ 
scription, the principle of a linear amplifier can be mod¬ 
eled by a second-tier of output voltage as v x = vo + xR , 
where R is the amplification coefficient and vq denotes a 
reference voltage. In addition to this one-to-one corre¬ 
spondence, the amplifier will inevitably introduce extra 
noise. Without loss of generality, we assume that the 
noise is Gaussian. Then, the output voltage v satisfies 
P x (v) = (27rtj 2 )~ 1 / 2 exp[—]. Viewing the relation 
between v x and x, we may reexpress this distribution 
alternatively as 


P x (v) 



(u — x) 2 
2a 2 


( 21 ) 


where we have re-defined a = a/R and v = (v — vq)/R. 

Now we analyze the consequence of the amplifier’s 
noise. Let us consider first how the initial probability of 
the first-tier output, Pi{x) = pnPi(x)+p 22 P 2 (x), should 
be modified from the amplified signal v (or the re-scaled 
v) . This can be done by using the classical Bayes formula: 


P v (x) 


Pi{x)P x {v) 
f dxPi{x)P x (v) 


( 22 ) 


This distribution tells us that, for a given amplified sig¬ 
nal v, there are many possible x. Therefore, the re¬ 
conditioned qubit state should be 


Pv — 


dxp x P v (x ) , 


(23) 


where p x is the qubit state conditioned on x, in general 
which is given by the quantum Bayesian rule of Eq. (js|). 
More explicitly, we have 

Pv,n = pi\Bi{v)/N v , (24a) 


Pv,22 — P22B>2{y)l Af v , (24b) 


Pv, 12 — PuB^{v)/M v ■ (24c) 


Here, for the sake of brevity, we have introduced 

(„-xi, 2 ) 2 (a _.2 (v-xnA 

Bi, 2 {v) = e 2 < d+ct2 > , Bz(v) = e~ sn e 2 ( D +^ 2 ) and 
A f v = pnBi(v) + P 22 B 2 (^); where Ax = \x\ — x 2 \ and 
xo = (Si +S 2 )/2. 

Similar as using Eq. © , we can calculate now the 
“noisy” WV by means of 


_ / dwPi(v)P v (f) 
f[X)i f dvP(v)P v (f) 


(25) 


where Pi(v) is the probability distribution of the second- 
tier result “v” which is given by 


Piiv) = J dxPi{x)P x {y) = 


PuBiiy) + P22B2 (u) 
yj2-K(a 2 +DR 2 ) 


(26) 


and the post-selection probability reads 

Pv(f) = Pfllpv,ll + Pf22pv,22 + 2 Re(p/I 2 p u , 2 l) • (27) 


In this way we obtain the result of WV defined from 
the “polluted” data, after noisy amplification. However, 
being of some surprise is that we find the same result 
of Eq. (|l 8 ]). This means that the WV is free from the 
amplifier’s noise. We may understand this result by the 
fact that the noisy amplification of the second-tier does 
affect the distribution width of the outputs, but does not 
affect their average even in the presence of post-selection. 

The extra noise introduced in the amplification process 
is a sort of non-ideality to reduce the quantum efficiency 
of measurement. One may expect and actually can prove 
that other sources of non-ideality owing to quantum in¬ 
formation loss does not affect the WV as well. Desirably, 
the feature that the WV is free from amplifier’s noise 
and other sources of non-ideality (quantum information 
loss) can benefit the measurement and applications of the 
WVs. 


III. WEAK VALUES IN CIRCUIT-QED 


In this section we specify the WV study carried out 
above to the solid-state circuit-QED (cQED) system for 
two reasons. First, the cQED system is one of the most 
promising solid-state architectures for quantum informa¬ 
tion processing [30, ^ and excellent platform for quan¬ 
tum measurement and control studies jl2] [m). Second, 
this is to-date the most experimentally accessible solid- 
state system where the continuous weak measurement of 
the type considered in Sec. II has been realized [|2] [34| . 
In particular, we will obtain even more useful result than 
Eq. © in that it allows for direct measurement of the 
full AAV WV, thus for efficient method of qubit state 
tomography in this important system. 
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A. QTE Approach 


Under reasonable approximations, the cQED system 
is well described by the Jaynes-Cummings Hamilto¬ 
nian fiof . Moreover, in dispersive regime, by perform¬ 
ing a qubit-state-dependent displacement transformation 
(called also “polaron” transformation), it is possible to 
eliminate the cavity degrees of freedom to get a trans¬ 
formed QTE for the qubit state alone under continuous 
homodyne measurements ||f4|: 


.£jq-j-B(t) 1 I\j(£) r 

p = - - [a- ,p] + ^-V[a }p 


-y/f ~(f)M[a z ]pZ(t) + ■ 


(28) 


Here, to q = uj q + y is the renormalized qubit energy owing 
to a dispersive shift % to the bare energy u> q ; B(t) = 
2yRe[ai(£)a|(£)] is a generalized dynamic ac-Stark shift, 
with ai(t) and a-i(t) the cavity fields associated with the 
qubit states |1) and |2), respectively. In addition to the 
Lindblad term, a new superoperator is introduced here 
by M[a z }p = (a z p+pa z )/2— (a z )p, with ( a z ) = Tr [cr z p\. 
Also, 


r d(Q = 2xIm[ai(£)Q2(£)] , 
r d{t) = n\/3(t)\ 2 cos 2 (ip - Op ), 

T ba (t) = K,\/3(t)\ 2 sm 2 (tp - Op), (29) 

characterize, respectively, the ensemble-average dephas¬ 
ing, information-gain, and back-action rates. In these 
expressions, we have introduced (3{t) = a. 2 (t) — a\{t) = 
\P(t)\e lS/3 , and denoted the cavity damping rate by k and 
the local oscillator (LO) phase by p for the homodyne 
measurement. 

Within the framework of “polaron” transformation, 
the homodyne current of the cavity-field-quadrature mea¬ 
surement is reduced to J(t) = — y/1\i{cr z ) + £(£). As ex¬ 
plained in Sec. II, from this result we can extract the 
Wiener increment via dW = x + y/T C i(a z )dt, where 
x = J(t)dt is the output of quadrature measurement. 
Accordingly, we calculate the post-selection restricted av¬ 
erage of “x” (the weak value) as follows. First, with 
the help of the above QTE, we update the qubit state 
from the “initial” state p(t), which is denoted also by 
Pi = t° p(t + dt) based on the measurement re¬ 

sult “x”. Then, we calculate the success probability of 
post-selection via P x (f) = {ipf\p(t + dt)\ipf) which yields 


P x (f) = P 1 + 


H + U)ImK) + ^(K| 2 -l) 


dt 


Vr“(Re(a^) - (a z )) + y/r^Im«,)] dw} , (30) 


and carrying out the integration of f dxxP(x)P x (f) and 
J dxP(x)P x (f), we obtain 

/v = _ ei Re(cr^) + e 2 Im(<7*) _ 

f l + [nim(a z w ) + (T d /2)(\a z w \^-l)]dt' 

where ei = y/T c i dt , £2 = y/T^dt and H = co q + B. In 
deriving this result we have identified Xy 2 = T V^cidt 
and D = dt, similar as in the previous section. 

Eq. ( |3l| ) is a useful result. The particular structure in 
the numerator, i.e., the presence of both the real and 
imaginary parts of allows for convenient measure¬ 
ment of the full complex AAV WV. Indeed, from the ex¬ 
pressions of T c j and Tf, a , we can selectively make either 
T ci or T;, 0 be zero, by tuning the LO’s phase ip. Then, 
in short-time limit (which in most cases makes the de¬ 
nominator near unity), one can separately measure the 
real and imaginary parts of the AAV WV (a^). A little 
bit complexity is for unknown initial state \ipi)- In this 
case, the post-selection of state | ipf) cannot rule out re¬ 
sulting in a large a ^, which would make the denominator 
of Eq. (|3l]) considerably deviate from unity. In practice, 
one can solve this problem by iteratively substituting the 
imprecise a z w into the denominator of Eq. ([nj) to obtain 
better estimation for Convergence is expected by a 
few such iterations. 

The ability of measuring the real and imaginary parts 
of the complex AAV’s WV has important applications. 
One of them is to develop new technique of state tomog¬ 
raphy ^-[h]]. For such purpose in cQED system, one 
should identify the ac-Stark shift and all the rates (T C j, 
Tf, a and T^) in Eq. (|dl]). From Eq. (^9|) we see that in 
general these quantities are of time dependence, i.e., de¬ 
pending on the temporal evolution of the cavity field. In 
experiments, however, the cQED setup is usually pre¬ 
pared in the “bad-cavity” and weak-response limits 
p~i| . [l8t . In this case, the cavity-field would evolve to sta¬ 
tionary state on timescale considerably shorter than the 
quadrature measurement time (denoted in this work by 
11 dt" ) . Therefore, one can simply obtain the ac-Stark 

shift and all the rates in Eq. ( |3l| ) by using the stationary 
coherent-state fields d\ and d 2 , which read 

di( 2 ) = -ie m /[-i{A r ± y) + k/ 2], (32) 

where A r = uj m — uj r is the offset of the measurement 
and cavity frequencies. For instance, in the bad-cavity 
and weak-response limit, we obtain the stationary B(t) 
as B ~ 2yh, which recovers the standard ac-Stark shift 
by noting that n = |d| 2 and a = —ie m /(^). Also, for 
resonant drive (w m = ia r ), we have 9p = 0. Therefore 
one can measure the real and imaginary parts of the AAV 
WV by choosing ip = 0 and 7 t/ 2, respectively. 


B. Bayesian Approach 


where P \(ijjf |^»)| and a w ’ as i n ^ 10 Let us continue to consider the WV for finite strength 

duced above. Further, replacing dW with x + y/T C i(a z )dt, measurement in cQED, applying the Bayesian scheme 
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proposed in Refs. (4|, [ifj. In general, for finite time 
( t m ) quadrature measurement, the output corresponds to 
x = x{tm) = So m dtJ(t ), where J(t) = -VT^(a z ) + £(i). 
Associated with qubit state |j), the priori knowledge of 
distribution Pj (x) of the integrated quadrature x is Gaus¬ 
sian, given by Pj{x) = -?== cxp[— (x — Xj) 2 /(2D)\ where 

Xj = ( 1 )'^ F cl tm and D = tm- 

With the knowledge of Pj(x), the quantum Bayesian 
rule is the same as Eq. (||) for the diagonal elements, how¬ 
ever for the off-diagonal element it needs the following 
essential corrections Jfffl 

Pl2(jm) ~ Pl2{tm) | (c^2 (tm) | (^m)) | 

x exp{-*[$i(f m ) + $ 2 (*m)]} • (33) 

In this expression, pi 2 (t m ) is the updated result given 
by Eq. (||). The overlap of the cavity-field coherent 
states, \{ot 2 {tm)\<xi(t m ))\, characterizes the effect of pu¬ 
rity degradation. Other corrections are two phase factors: 

$i(*m)= f B(t) dt , (34a) 

Jo 

$2(tm) = - f m V^ba(t) J(t) dt . (34b) 

Jo 

The first factor stems from the energy shift of the qubit 
owing to dynamic ac-Stark effect, and the second one 
is from the last term of Eq. (^) . A tricky issue in¬ 
volved here may need to be reminded. If we only con¬ 
sider the last term of Eq. (HI), the correction should be 
d >2 (t m ) = — / Q tm s/r ba (t) £(i) dt. However, for each step 
of evolution over (t,t + dt), the stochastic output current 
is J{t) which has a mean of J(t) = —^T ci (t)(a z )t- To 
guarantee the ensemble average of the off-diagonal ele¬ 
ment over the stochastic x = J(t)dt. to be valid within 
the Bayesian scheme, one should adopt J(t) rather than 
£(i) in the integrand of $2 (tm). This is consistent with 
Eq. (|2^) after ensemble-average over £(t)dt = dW(t), in 
that the last term vanishes and does not leave any phase 
factor such as e - l V / GW)A t )< it _ 

Similar to the short-time measurement analysis, be¬ 
low we restrict our consideration in the “bad”-cavity and 
weak-response limit. In this case, the ac-Stark field and 
all the rates are independent of time, being given by the 
steady-state cavity fields of Eq. (|||) . As a consequence, 
the both phase factors are simplified to $1 (t m ) = Bt m 
and $ 2 (t m ) = — \/r c i x(t m ), and the purity-degradation 
factor \(a 2 (t m )\a.i(tm))\ can be approximated to unity. 
After these, the WV can be straightforwardly calculated 
using Eq. 0 ) , where the post-selection probability P x (/) 
depends on the updated state p(t m ) from the cQED 
Bayesian rule, particularly with Eq. (|3S|). 

Like Eq. (|l^), for cQED we obtain j(x)i = M\/M 2 , 
with Mi and M 2 given by 

Ml = ~\/r C i tm(PfllPll — P/22P22) 

+ 2^r^t m e~ rdtm Im(p /12 pl 2 e“ m ), (35a) 


M 2 = {pfiipn + P/22P22) 

+ 2 e~ Tdtm B.e(pfi 2 p* 12 e i ^ ltm ), (35b) 


where T^ = (T ci + T ba )/2 is the overall measurement- 
induced decoherence rate. From this result, one can eas¬ 
ily recover the WV of Eq. ([il]) for short-time limit. How¬ 
ever, even for general finite time measurement, we can 
obtain similar elegant expression as well. After some al¬ 
gebra based on Eq. (||i|), we get 


61 Re(g* ) + e 2 Im(d* ) 

1 + S(\di,\ 2 ~ 1 ) 


(36) 


where ei = y/T^t m , £2 = i/Pb^tme and Q = 

(1 — e~ rdtm )/2. Notably, <r* in this result is a slightly 


modified AAV WV, taking a form as d z = 

where |i pi) differs from the initial state \ipi) = ci 11)+C 2 12) 
by a phase factor in terms of |i pi) = cie _ifMm |l) + C 2 12). 

As a final remark, similar to the weak measurement 
limit discussed above in previous subsection, by tuning 
the LO phase <p based on Eq. ([ib]), one can still conve¬ 
niently measure the real and imaginary parts of o/j,, from 
which efficient state-tomography technique can be devel¬ 
oped for finite strength measurement in cQED system. 
In practice this will make the experiment easier since it 
allows for larger output signal. 


IV. CONCLUSIONS 

To summarize, by applying both the quantum- 
trajectory-equation (QTE) and quantum Bayesian ap¬ 
proach, we obtained non-perturbative and exact weak 
values (WVs) for continuous weak measurement of 
qubits. Differing from the usual unitary interaction 
model for weak measurement in WV studies, the con¬ 
tinuous weak measurement scheme allows for non- 
perturbative Bayesian approach to update the qubit 
state, yielding thus exact WV result. From this, by 
reducing the measurement strength (time duration), we 
obtained exact expression for WVs in the weak measure¬ 
ment limit, which is in full agreement with that derived 
by using the QTE scheme. 

In particular, we have also extended the study to 
circuit-QED system. The obtained desirable results al¬ 
low for convenient measurement of the real and imag¬ 
inary parts of the AAV WV, and thus for developing 
efficient technique of state-tomography. Note that the 
cQED is to-date the most experimentally accessible solid- 
state system where the continuous weak measurement 
considered in this work has been realized. Moreover, as 
analyzed in Sec. II B 3, the WV measurement is free from 
the amplifier’s noise and/or the quantum efficiency of 
measurement. Therefore, the cQED system is expected 
to be an ideal platform for WV studies and related ap¬ 
plications, such as developing technique of direct state- 
tomography. 
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